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ABSTRACT 

Recent studies of the tight scahng relations between the masses of supermassive black holes and 
their host galaxies have suggested that in the past black holes constituted a larger fraction of their host 
galaxies' mass. However, these arguments are limited by selection effects and difficulties in determining 
robust host galaxy masses at high redshifts. Here we report the first results of a new, complementary 
diagnostic route: we directly determine a dynamical host galaxy mass for the z = 1.3 luminous quasar 
J090543. 56+043347.3 through high-spatial-resolution (0.47", 4kpc FWHM) observations of the host 
galaxy gas kinematics over 30x40 kpc using ESO/VLT/SINFONI with LGS/AO. Combining our result 
of Mdyn = 2.05lJ;^^ X 10" Mq (within a radius 5.25 ±1.05 kpc) with Mbh.MsII = 9. 02 ±1.43 x 10* Mq, 
Mbh.Hq = 2.83li['^3 x 10* Mq, we find that the ratio of black hole mass to host galaxy dynamical mass 
for J090543. 56+043347. 3 matches the present-day relation for Mbh vs. MBuigG,Dyn7 well within the IR 
scatter, deviating at most a factor of two from the mean. J090543. 56+043347. 3 displays clear signs of 
an ongoing tidal interaction and of spatially extended star formation at a rate of 50 — lOOMoyr"^, 
above the cosmic average for a galaxy of this mass and redshift. We argue that its subsequent evolution 
may move J090543. 56+043347. 3 even closer to the z = relation for Mbh vs. MBuigc,Dyn- Our results 
support the picture where any substantive evolution in these relations must occur prior to z ~ 1.3. 
Having demonstrated the power of this modelling approach we are currently analyzing similar data 
on seven further objects to better constrain such evolution. 

Subject headings: galaxies: active - galaxies: kinematics and dynamics - quasars: emission lines - 
quasars: individual (SDSS J090543.56+043347.3) 
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1. INTRODUCTION 

The tight correlations between the properties of 
supermassive black holes (BHs) and their host galax- 
ies jK ormcnd v & Richstonc 1995' Magorrian et al. 
19981: ^ TCebhardt et all 120001: iFerrarese'fc MerritI 



200a iTremaine et all I2002t iMcLure fc Dunlod 120021: 



Marconi k HuntI 120031: iHaring fc Rixl 120041: iGraham 



20041: iKormend v fc Benderl l2009t iGiiltekin et al.' '200£ 
Jahnke et a'Ln 2009: Merlon i et all [2010: Graham ct alj 



20111) are a powerful tool for probing their relative and 
possibly interdependent evolution and growth across cos- 
mic time. Since the hierarchical assembly in conjunction 
with the central limi t theorem can explain the existence 
of these correlatio i is (lPeng|l2007l : iHirschmann et al]l2010l : 
IJahnke fc Macciol 1201 Ih . evolution in the correlations 
may solely be due to the relative growth speeds of 
the black hole mass and the corresponding total and 
bulge stellar masses of the host galaxy as a function 
of redshift. Hence, evolution in these correlations is 
a key diagnostic of the growth channels and fuelling 
mechanisms of black holes, active galactic nuclei (AGN) 
and their host galaxies. 

Overall, the current picture is one of evolution in 
th e various correlatio i is clearly being s tro ng at z > 
2 (l^yalter efall [2005 iPeng et al.l [2006allbl: iSomerville I 
120091 iDecarT et al T l2om A~ number of studies also 
find a strong evolution signa l at low redshif ts (jWoo et al.l 
[1]06; Trcu ct al. 2007; Ben nert et alll2010[) . particularly 
in terms of correlations involving the bulge properties 

|inskip@mpia.de] 



of the host galaxies: it is generally recognised that black 
holes of a given mass are hosted by increasingly less mas- 
sive spheroids/bulges at larger values of z. However, at 
lower redshifts the majority of studies show a mild (by 
factor ^2-3) or absent evolution, particularly where the 
total stellar mass of the host gala xy is concerned. For 
example, for AGN at 1 < z < 2 iMerloni et all (|20100 
observe an offset of a factor of 2-3 from the z = scaling 
relation, and note that given the large scatter, the data 
could still be co nsistent with zero e volution at the 2a 
level. And while IPecarli et aD ()2010t ) observe an overall 
evolution of up to a factor of 7 between z = and z = 3, 
the sources observed at intermediate (1.0 < z < 1.5) and 
low (z < 1) redshifts are consistent with modest (0.3dex) 
and no evolution respectively. 

Beyond the nearby Universe, observational estimates 
of black hole and galaxy masses are subject to consider- 
able uncertainties. Furthermore, different measures of 
galaxy mass (i.e. bulge mass, stellar mass, total mass) 
are frequently used, necessitating additional care in the 
interpretation of results from different studies. The 
major uncertainties of previous work in this field affect 
both large and small samples equally, and arise from the 
accuracy with which the masses of black hole and host 
galaxy (total and/or bulge mass) can be determined. 
Outside the local universe, black hole masses can only 
realistically be obtained for currently accreting type-1 
AGN (i.e. quasars or Seyfert-1 galaxies). For such 
objects, Mbh can be constrained from the combination 
of the UV-luminosity and the broad emission line 
widths from ionization models which assume virialised 
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motion and are in turn calibrated via reverberation 
mapping on local AGN with a statistica l uncertainty 
withi n a factor of 3-4 for Mgli a nd Civ (iVestergaard I 
^nol iMcLure fc .Tarvisl [2001 IShields et al.| (20?)^ 
Greene & Ho 2005; Vestergaard & Peterson 2006 
Woo et al . 2006; Kollnieier et al. .2006 ; Salyiander et ^alJ 
20071: iNetz er & Trakht ejArotl llool iTreu et al.l 126077 
McGill et al.li200&.Onken fc K ollmeier II 2008D . orl ess for 
Balmer lines pennev et al.l 12009; Shen fc Kellvl [20Tol ) . 



plus any further systematic uncertainty associated with 
the calibration between reverberation mapping BH 
masses and dynamical BH masses at z—0. 

However, obtaining accurate (total and/or bulge) dy- 
namical masses for the quasar host galaxy, the second 
key ingredient of the correlation in question, is very dif- 
ficult for objects beyond z > 0.5. The presence of an 
active nucleus makes stellar velocity dispersion measure- 
ments from absorption lines infeasible, and the estimated 
velocity dispersions derived from narrow emission line 
widths (IShields et al.l[2003l[2 00l 'Salviander et al.l[2007l: 
iHo et al.l I2008t iRiechersgr a l. 2008, 2003) are f raught 
with uncertaintv ()Ho 1 120071 : iGreene et all 120091 ). In- 
ferring stellar masses from imaging-based measurements 
of the host galaxy stellar luminosity is another pos- 
sibl e approach use d by a number of pre vious stud- 
ies (iBorvs et al.|[200 5: Pen g et al.|[206 6a.b; J ahnke et ahl 
120091 iSomerville 1 120091; iMerloni et al.li201Q) but, as the 
mass-to- light ratios deduced from the UV/optical spec- 
tral energy distributions of star-forming galaxies are 
severely model-dependent, inferred stellar masses are 
non- unique/uncertain. At redshifts z > 4, practical 
mass determination is limited to methods such as the 
dete rmination of dynam ical masses via CO observations 
(e.g. [Walter et al.l 1200^ . Direct measurement of a host 
galaxy's dynamical mass is free of many of the problems 
associated with other methods, and can equally be ap- 
plied to observations at lower redshifts. Aside from mm- 
observations of CO emission, dynamical masses can also 
be measured from ISM emission lines in the rest-frame 
optical. 

At redshifts of z > 1 , Ha is the most prevalent power- 
ful emission line that can be used in practice as a kine- 
matic tracer in AGN host galaxies. Even so, successfully 
probing the geometry and kinematics of the Ha-emitting 
gas requires integral field near-infrared (near-IR) spectro- 
scopic observations with a spatial resolution of not more 
than a few kpc (i.e. < 0.5"). A sufficiently high (nar- 
row line) luminosity is needed, but the Ha luminosity of 
galaxies undergoing star formation at the average rate is 
easily sufficient - along with the presence of favourably 
orientated and ordered velocity structures (also likely for 
a significant proportion of objects). As a pilot study for 
direct dynamical host galaxy masses for z > 1 quasars 
we present and model near-IR integral field spectroscopy 
of the z =1.3118 quasar J090543.56-F043347.3 using SIN- 
FONI ([Bonnet et al.|[20p; lEisenhauer et al.l 120031 ). We 
intend to present the results of similar observations of 
our full sample in a subsequent paper; interim analysis 
suggests that the required fiux, orientation and ordered 
velocity field constraints are met for approximately 2/3 
of the objects targeted. A standard cosmological model 
with n = 0.27, A = 0.73 and Hq ^ 71kms-iMpc-i 
is assumed throughout; this leads to a spatial scale of 
8.45 kpc/arcsec at this redshift. 



2. OBSERVATIONS AND DATA REDUCTION 

Our target, the z = 1.31177 SPSS quasar 
J090543.56-f'043347.3 ([Schneider et al.l[2007[ : [Shen et al.l 
[200 81. was selected on the basis of the presence of a 
nearby guide star (12th magnitude at a distance ^^4.3 
arcmin) to be used for tip-tilt correction. This ob- 
ject was observed using the near-infrared integral field 
spectrometer SINFONI on the Very Large Telescope of 
the European Southerii Observatory (Bonnet_et al. 20O3; 
[Eisenhauer eFall [20Q1 and the PARSEC Laser Guide 
Star Adaptive Optics, on the night of 2009 December 
15, with an average airmass of ^1.16 and IFU image- 
seeing of 0.47arcsec FWHM. Six pointings of 600s were 
obtained, in two repeating patterns of target-sky-target 
in order to facilitate the removal of sky background emis- 
sion, with a sky-frame declination offset of 30 arcsec from 
the target coordinates. The H-band filter was used with 
a grating central wavelength of 1650nm, giving a spec- 
tral coverage of 1.45-1.83/im and a spectral resolution 
of 3000. The spatial scale of each of the 32 SINFONI 
slitlets was lOOmas; these are imaged onto 64 detector 
pixels apiece, giving a field of view of 3" x 3" for each 
pointing. 

Analysis of SINFONI data requires a painstaking treat- 
ment of cosmic ray hits and bad pixels, neither of which 
can be adequately managed by the existing pipeline re- 
duction routines alone. We use IRAF routines to identify 
bad pixels in each raw frame of data prior to applying the 
instrument pipeline, creating a mask file for each individ- 
ual raw frame. Otherwise, the initial ^2-3% bad pixels 
fiagged as NaN in each raw frame would breed dispro- 
portionately as the raw data are straightened, stretched, 
shifted and combined within the pipeline, resulting in an 
unrealistic fraction of the final datacube containing NaN 
values (typically '--^30-50%). Ideally, we do not seek to 
remove signal from any pixel with some level of contam- 
ination by bad data, but rather to trace the propagation 
of bad pixels and thus gauge the reliability of each data 
element within the final datacube. From this point on- 
wards we work with two different data sets. Firstly, we 
use our final bad pixel mask to linearly interpolate over 
all flagged bad pixels in the observed data. Although 
these flagged pixels will contribute flux to their neigh- 
bours, they will only do so with locally realistic count 
values. Hence, we do not invent data while preserving the 
integrity of the spaxel grid. Secondly, mask frame images 
are generated and treated in the exact same manner as 
the observed data, thus allowing us to keep track of the 
behaviour and growth of the bad pixels as the datacube 
is generated. 

A further issue is the mis-location of slitlet edges by the 
ESO SINFONI data reduction pipeline. To counter this, 
we apply small pixel shifts to the supplied WAVE, MAP 
calibration images, recovering 2-3 of the 3-4 columns of 
data on each side of the held of view i.e. 10% of the held, 
which would otherwise be lost. The outermost columns 
of each slitlet which cannot be reliably wavelength cali- 
brated remain masked out. Each raw frame and its as- 
sociated mask are then individually converted from a 2-d 
image into 3-d datacubes using the 'sinfo-rec-jitter' rou- 
tine in the pipeline and the improved skyline subtraction 
IDL routines of Davies (2007). Our bad pixel mask dat- 
acubes contain the proportional contribution of bad pix- 
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Fig. 1. — From datacube to velocity field, (a - left) Illustration of the continuum-subtracted QSO spectrum at the peak spaxel. The 
two wavelength regions of the broad Ha line found to be optimal for our QSO PSF construction are highlighted in grey. They are free of 
contamination from narrow Ha or [Nil] line emission, whether it be due to the AGN itself or the spatially extended emission from the host 
galaxy, or the night OH sky lines. While there are 3-4 other possible wavelength regions within the broad Ha line which are uncontaminated 
by narrow line emission or night sky lines, the flux in these regions has lower signal-to-noise, and results in an overall worsening in quality 
of the resulting PSF profile, (b - right) The resulting PSF image defined as the spaxel by spaxel intensity map in these three wavelength 
regions. Contours are at 1, 2, 5, 10, 20 and 50% of the peak QSO flux, with the lowest (1%) contour in black. The 50% and 70% encircled 
energy radii are 0.12" and 0.2" respectively, corresponding to contours at 40% and 14% of the peak flux. 



els to each spaxel in the cube. For spaxels where the pro- 
portional contribution of bad pixels is larger than 50%, 
that particular spaxel in the corresponding observational 
data cube is converted to a value of NaN; such spaxels 
do not contain reliable information. This accounts for a 
total of approximately 6% of spaxels in our final com- 
bined datacube. The remaining spaxels in the datacube 
are either completely uncontaminated by the presence of 
nearby bad pixels (true for ^^50-70% of spaxels) or are 
dominated by emission from clean spaxels (~25-45%). 
As described earlier, in this latter case the spaxel flux is 
derived from an initial image where the flagged bad pix- 
els/cosmic ray hits have been replaced by linear interpo- 
lation over the surrounding clean data. As these spax- 
els originate from a majority of good data pixels, with 
a minority contribution from proportionately scaled lo- 
cally plausible flux values, they are judged reliable, and 
left as they are. This done, the individual datacubes 
are straightened and aligned via a process of gaussian 
centroiding, sub-pixel linear interpolation and shifting of 
each wavelength slice. 

Our final datacube is produced by mean-combining 
the individual aligned cubes on a spaxel-by-spaxel ba- 
sis, ignoring masked NaN values and 3-sigma clipping 
any value which deviates significantly from the typi- 
cal values found at that spaxel position and those of 
the neighbouring spaxels across all individual aligned 
cubes. Using SINFONI observations of the standard 
stars Hip047903 (GIV) and Hip047235 (G3V), and tem- 
plate spectra (Pickles 1998) for their specific spectral 
types scaled to the appropriate magnitudes in the H- 
band, we determine the relative flux response for each 
wavelength slice across the datacube, and apply the de- 
rived correction (accurate to ~ 12% in absolute flux val- 
ues) to the cube data. 

3. ANALYSIS AND RESULTS 



3.1. The spatially resolved narrow Ha emission 

In order to reveal the underlying kinematics of the host 
galaxy, we need to extract the extended narrow line emis- 
sion from our datacube. We firstly remove the contin- 
uum emission in each spaxel via a least-squares fit to a 
second-order polynomial (excluding the wavelength re- 
gion covered by the broad Ha emission), and then turn 
to the more difficult task of separating the line emission 
of the host galaxy from that of the AGN itself. As the 
broad lines originate only from the (unresolved) imme- 
diate vicinity of the AGN, the observed spatial variation 
of broad line intensity in our datacube can be used to 
accurately constrain the point spread functi on for our 
data (as demonstrated bv lJahnke et al.ll2004[ ). For each 
spaxel in our data cube, we isolate the wavelength ranges 
of the spectrum which contain strong broad line emis- 
sion, but without contamination by sky line residuals, 
to either side of the spatially extended but spectrally 
compact narrow line emission (Fig.[T^). Using the emis- 
sion either side of broad Ha, we subtract any residual 
background continuum level, and sum up the remaining 
flux in the selected regions of the broad line. This pro- 
vides us with a relative measure of the broad line flux 
on a spaxel-by-spaxel basis in the form of a PSF image 
which accurately traces the relative distribution of light 
from the quasar point source at all wavelengths within 
our datacube (Fig. ^p). The quasar spectrum itself is 
well described by the spectrum of the peak spaxel in our 
datacube, which is dominated by AGN emission alone. 
We construct a quasar-emission datacube by scaling the 
quasar spectrum by the relative flux in each pixel of our 
PSF image. Subtraction of this datacube from our cal- 
ibrated data cleanly and accurately removes both the 
broad and narrow Ha emission due to the AGN, and re- 
sults in a new datacube which contains only host galaxy 
emission. Non-negativity within the wavelength region 
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Fig. 2. — Example fits to the narrow fine emission in the post- 
QSO subtraction datacube for tiiree different spaxels, one in tfie 
tidal feature to the south east (top), one to the north (center), 
and one in the central region (bottom). The labelled positions are 
in the same coordinates as shown on Fig. lb and all subsequent 
2-d images of the field of view. Gaussian line profiles plus a lin- 
ear continuum are used to simultaneously fit the narrow Ha and 
[Nll]6548,6584 emission lines in the source rest-frame, assuming the 
canonical SDSS DR3 redshift of 2 = 1.31177. The weaker [Nil] lines 
cannot be accurately identified in all spaxels, and in such cases are 
excluded from the fit. 

covered by the broad Ha emission is not enforced. For 
twenty spaxels close to the center of the PSF, errors in 
the subtraction of the highly luminous QSO emission 
dominate the resulting spectra and in a few cases result 
in mild over-subtraction of the QSO emission. We do not 
retrieve reliable information on the narrow Ha emission 
from these spaxels. 

For each spaxel of the QSO-subtracted datacube, we 
then fit the data with a linear continuum and a simulta- 
neous combination of three gaussians of the same width 
and fixed relative centroid but varying strength for the 
Ha and [Nii]6548A,6584A emission lines. The spectra 
are converted to the source rest frame using the canon- 
ical SDSS DR3 redshift of z = 1.31177. Fig. [2] displays 
example fits to three different spaxel spectra in the QSO- 
subtracted datacube from different regions of the field of 
view. The results are then masked to exclude unreliable 
fits to noise rather than to genuine line emission, i.e. ex- 
cluding the outermost rows and columns of data where no 
genuine emission features are observed, spaxels where the 



fit results in extreme velocities (—100 < v < 600km s~^ 
without subtraction of any further systemic velocity at 
the assumed redshift of 1.31177, confirmed as fits to noise 
features via visual inspection), line centroid uncertainties 
> 40km s~^, and Ha fluxes of < 10 counts or uncer- 
tainties of > 3 counts per binned spaxel. Fig. [3] dis- 
plays the extracted Ha flux image (Fig. [5^), line width 
image (Fig. [Sb), velocity field after correcting the sys- 
temic velocity (Fig. HJ;) and [Nii]6584A/Ha fine ratio 
map (Fig. [3ji), overlaid with the PSF contours from 
Fig. [TJj. A more in-depth analysis of the emission line 
properties of this system is deferred to our study of a 
larger sample of similar objects. In our subsequent mod- 
elling of the velocity field, we further mask the outermost, 
low signal-to-noise regions, via boxcar-averaging the line 
fiux in 9 by 9 pixel regions and excluding regions of low 
average signal. The remaining unmasked pixels are illus- 
trated in the line fiux image displayed in Fig. 2^. 

3.2. Velocity field analysis 

As is clear from Fig. [5};, this source displays a com- 
plex velocity field. A quick visual inspection suggests the 
central region can be well described by ordered disk-like 
rotation of the gas, with a reversal of the gas velocity at 
larger radii along at least one arm-like feature, and possi- 
bly a second outside our field of view to the north. With 
a pixel scale of ^ 0.4 kpc per pixel, the southern arm 
extends to a projected distance of approximately 20 kpc, 
suggesting a tidal origin for this feature as the most likely 
scenario. Although there are no radio observations made 
specifically of this object, the fact that it is not a FI RST 
([Becker. White & Helfand 199^; iWhite et all Hw) ra- 
dio source suggests that no strong radio jet is present in 
this source, and thus that jet-driven shocks/interactions 
are unlikely to be a significant factor for the observed 
kinematics of this object 

Our aim is to use this observed gas velocity field as a 
tracer of the overall gravitational potential of the galaxy. 
The fact that disk-like rotation is present in the gas ve- 
locity field, as we will show in more detail in the fol- 
lowing sections, is independent of the geometrical struc- 
ture of the underlying stellar body and makes no state- 
ments about the stars being in either a disk-dominated 
or bulge-dominated structur^B- Gas di sks can be com- 
monly obse rved in elliptical g alaxies (jOosterloo et al.l 
I20ia 120071: lEmonts et alll2008l) . and it is often the case 
that their kinematics do not trace those of the under- 
lying stellar population. Indeed, the velocity structures 
may be highly decoupled, both in the case of active (e.g. 
HE1029-1408 at z = 0.086, which dis plays low stellar 
rotat ion but substantial gas rotation; iHusemann et ahl 
I20T0I) or inactiv e (e.g. NGC2768 NGC4314, NGC 4526 
and NGC 4546; ISarzi et al.l[2006l) galaxies. 

In order to properly characterise the rotation of the 
central regions and to extract a dynamical mass, it is 
important that we can also simultaneously account for 
the contrasting peculiar velocity structures dominating 
the gas velocity field at larger radii (> 5kpc). However, 
the signal-to-noise of the line emission in the individual 
spaxels where these two velocity structures overlap is not 

^ We note that in spite of the depth of our observations, there is 
insufficient extended continuum flux within the datacube to extract 
any information about the host galaxy morphology. 
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Fig. 3. — Results of emission line modelling (see section 3.1), overlaid with PSF contours from Fig. lb. (a - top left) Narrow Ha emission 
line flux image, (b - top right) Line width of modelled narrow line emission, (c - bottom left) Narrow line emission velocity field. Note 
the change in the rotation direction between the inner and outer regions. Zero velocity is defined by the center of rotation in the central 
regions, (d - bottom right) [Nll]6584/Ha emission line ratio. 



sufSciently high to separately resolve more than one line 
profile, nor can multi-spaxel extractions allow the fitting 
of multiple line profiles without unavoidable degeneracies 
in the relative line widths, centroids and strengths. We 
therefore choose to use multi-component fits to the full 
velocity field. 

3.2.1. Separation of emission from the central and outer 
tidal components 

In the process of producing a model velocity field, it is 
essential to consider the strength of the contribution of 
line flux at different velocities from adjacent pixels during 
convolution with the observational PSF. The observed 
velocity in a given position is not simply the convolution 
of the model velocity field with the observational PSF 
centered on that location: as the line flux varies across 
the observational field of view, the proportional contribu- 
tion of different pixels to each specific value in the output 
model velocity field will also vary. In order to correctly 
weight each velocity component of our model during PSF 
convolution, we also need to include a model of the rela- 
tive pixel-to-pixel line flux - for both the central rotating 
disk and also the additional tidal arm features. 

In producing a model for the velocity fleld of 
J090543. 56-1-043347. 3, we start by examining the two 



major velocity structures separately, followed by more 
extensive multiple-component modelling. We ob- 
tain additional constraints on these velocity structures 
via the ap plication of harmonic decomposition (cf. 
Ivan de Ven & Fathi 2010; Fathi et al. 2005). This treat- 
ment divides a velocity field and associated signal-to- 
noise map into elliptical annuli of differing semi-major 
axes but identical orientations, inclinations, ellipticities 
and centroids, and iteratively extracts the velocity pro- 
file as it appears within each annulus as a finite number 
of harmonic terms. Velocity and higher-order moments 
can then be readily characterised, and the nature of any 
deviations from non-circular motion assessed. The same 
tools can also be used to model a velocity field within 
a succession of elliptical annuli, while also varying the 
parameters listed above for each annulus. 

We use this technique to confirm the centroid of our 
velocity field, and find that the central regions are indeed 
well matched by ordered rotation, oriented at an angle 
of ~ 10° east of north, and inclined to the line of sight 
at an angle of 35 ± 15°. The results of fitting the ve- 
locity field with a single simple inclined rotating disk, of 
fixed position angle and inclination angle, and a veloc- 
ity which varies with radius, are displayed in Figs. |3^-d. 
These plots illustrate the velocity field, model, residu- 
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Fig. 4. — Projected velocity field and empirical first-order modelling. Here we examine separately the two major velocity structures in 
the emission line gas, prior to more extensive multiple-component modelling (Figs. 5-7). Harmonic decomposition techniques are used to 
fit circular rotation along two sets of ellipses of fixed position and inclination angles (see section 3.2.1). We display the data (frame a), the 
overall model velocity field for the first set of simple tilted ring models for which the position angle was tuned to target the central rotating 
region (frame b), model-subtracted residuals (frame c) and the resulting velocity as a function of radius for each extracted elliptical annulus, 
for this first set of fits (frame d). Rotational velocity increases to a maximum of 220 ± 50kms~^; at larger radii (beyond a radius of ~ 12 
pixels) it then declines steadily due to the infiuence of the counter-motions in the tidal features. Panels in the bottom row display the 
data (frame e), the model velocity field for the second set of fits with a position angle tuned to target the tidal feature (frame f), residuals 
(frame g) and the model velocity profile as a function of radius (frame h). In the region dominated by the tidal features, velocity varies 
roughly linearly with radius. 



als, and radial rotation profile of the model respectively 
(note that each data point in Fig. |4ji is derived from all 
spaxels which contribute to the model elliptical annulus 
at that radius). Velocity increases rapidly in the cen- 
tral regions to a maximum of 220 ± 50km s~^, and then 
declines rapidly beyond a radius of ~ 12 pixels due to 
the influence of the counter-motions in the tidal features 
(Fig.ffll). 

We also carry out similar simplistic modelling of the 
velocity field targeting the outer regions. The velocity 
field due to the tidal arm features is oriented at an an- 
gle of approximately 140° offset from the central struc- 
tures, varies relatively monotonically with radius, and 
dominates the velocity field beyond a projected radius 
of approximately 5kpc (= 12 pixels; Fig. H^-H^). Full 
harmonic decomposition along each set of the defined 
ellipses implies that neither the central rotating struc- 
ture nor the outer tidal arms show any strong evidence 
for more complex motions within the regions they re- 
spectively dominate; the inclusion of any further velocity 
structures in our multiple-component modelling is there- 
fore not required. 

While harmonic decomposition provides a rapid means 
of investigating the different types of velocity structure 
contained with the velocity field of this galaxy, and de- 
termining better constraints on the system's inclination 
angle, we still need to model both structures simultane- 
ously, particularly in order to determine the underlying 
rotational velocities in the overlap region. We model the 
light distribution of the tidal arm features by first tracing 
the observed line fiux with a logarithmic spiral (Fig. [5^). 
The average cross-sectional flux profile is determined and 
applied (Fig. Od), and subsequently convolved with the 
PSF (Fig. [TJj) to create a good match to the observed 



fiux (Fig. [5t) . We then model the residual fiux (Fig. [SJi) 
with an inclined exponential disk (oriented to match the 
central ordered rotation region of the velocity field) , also 
PSF-convolved (Fig. [SJb), to provide a first approxima- 
tion to the fiux-profile of the rotating disk. (Note that 
the residual fiux after tidal-arm subtraction was better 
described as an exponential disk than as an n > 2 Sersic 
profile, as expected for a rotating gas structure.) Fig. [ST 
displays our combined relative flux model for the disk and 
tidal arm features, unconvolved with the observational 
PSF. While this does not perfectly match the observed 
flux proflle of the narrow line emission, it does provide 
a very good approximation to the relative pixel-to-pixel 
contributions to the velocity field. 

3.2.2. Velocity field model and dynamical mass 

The velocity profile of the arm features is modelled 
as a linear change in velocity with location along the 
arms, defined from the outer regions free of any influence 
from the central rotating disk, and extrapolated inwards 
from there. For the central rotating disk, we assume for 
our model velocity field a Persic universal rotation curve 
(iPersic et al.l 119961 ) as a generally valid example of the 
rotation curve for a disk. This can be described as 



Vix) = V{R. 



1.97x1-22 




-1.6 exp[-0.4(L/U)] 



where x = R/Ropt- We consider models with 0.1 < 
L/Li, < 5.0, 'optical' radii Ropt (roughly three times the 
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Fig. 5. — Modelling the Ha flux distribution (see section 3.2.1) as a prerequisite to full flux-weighted and PSF-convolved multiple- 
component velocity modelling (results of which are displayed in Figs. 6-7). (a - top left) Narrow Ha emission line flux image, plus 
logarthmic-spiral flt tracing the tidal arms, (b - top center) Model cross-sectional flux proflle for tidal arms, (c - top right) Convolution of 
cross-sectional flux proflle with observational PSF, providing best-fit model to narrow line flux in tidal arms, (d - bottom left) Residual 
narrow H-alpha emission line flux image, after removal of modelled flux in the tidal arm features, (e - bottom center) Convolution of 
n=l sersic disk model with the observational PSF, within the orientation and inclination constraints supplied by the velocity fleld, which 
provides the best-fit model to the residual line fiux. (f - bottom right) Combination of disk and tidal arm flux models, unconvolved with 
observational PSF, used in flux weighting of velocity models. 

models as defined by the parameter bounds above. De- 
rived from these curves using the gas motions as a tracer 
of the overall potential, we also plot the enclosed dynam- 
ical mass as a function of radius (Fig. This latter 
plot assumes that Vrot = Vdrc- In reality, non-circular 
gravitational motions of the gas will also be present (e.g. 
streaming and turbulence in the ISM), and the derived 
dynamical mass profile is, strictly speaking, a lower limit. 
However, as is clear from Fig.tSb, where the line emission 
originates solely from the centrally rotating regions - the 
eastern and western quadrants of the centrally rotating 
structure, which are not subject to possible contamina- 
tion by the tidal features - the observed line widths are 
at their lowest and in fact are barely resolved above the 
instrumental resolution (~ 70kms~^). The intrinsic line 
dispersion due to non-circular gravitational motion in the 
gas will be minimal, and thus Vrot ~ Vdrc is a good ap- 
proximation. In Fig. [3 we display the masked data, best 
fit model and model-subtracted residuals for the velocity 
field. 

We can reliably trace the galaxy rotation curve out 
to a distance of between 10-15 pixels, equivalent to a 
galactic radius of 4.2 - 6.3 kpc. (Note that the 70% en- 
circled energy radius for the point spread function of our 
datacube is ~ 0.2", equivalent to a projected distance 
of ^ 1.7kpc at the redshift of this source.) Ideally, we 
would determine a dynamical mass with reference to the 
scale length of the stellar light in the host galaxy, but no 
such information is available for this object. (While we 



disk scale length, and defined by iPersic et al.l ()1996[ ) to 
include ~ 83% of the light from an exponential disk; in 
practise our models trace the velocity field out to twice 
the optical radius) ranging from 0.11-15 pixels (with a 
plate scale of ^0.42kpc/pixel at this redshift). We also 
allow the inclination and orientation of the rotation axis 
to vary within our previously determined constraints 
from harmonic decomposition and from a prior model 
using a rotating disk alone without simultaneously ac- 
counting for the tidal arm structures. For each pixel, the 
contribution from the two model components within a 21 
by 21 pixel region equivalent to the size of the observa- 
tional PSF are weighted according to their relative fiuxes 
and convolved with the PSF. 

Our best fit model is determined using a Levenberg- 
Marquardt least-squares minimization technique 
(|Markwardtl[2009l) . and has L = 0.25L^, Ropt = 10.7 



4-3.0 
-1.6 



and an inclination angle of 47 ± 12°. While the formal 
la parameter errors for our best fit model are 0.2 pixels 
for Ropt and 1.1° for the inclination angle, we note 
that reasonable fits to the data can also be obtained 
with inclination angles ranging from 35° to 59° before 
gross mis-fitting of the data becomes clear in the 
model-subtracted residuals. As the assumed inclination 
angle has the greatest impact on the derived dynamical 
mass, the full parameter errors given above are derived 
after applying these limiting values for the inclination 
angle, and account for the loosest possible constraints 
on disk inclination for this source. 

Fig. [6K displays the de-projected rotation curve of the 
central gas disk component as a function of radius in the 
rest frame for our best-fit model. We also display depro- 
jected rotation curves for our upper and lower limit mass 



have determined a best fit measure of Ropt ^ 10.7 pixels 
for the gaseous emission for this galaxy, this does not 
necessarily have any bearing on the scale length of the 
stellar light, and should not be used as a proxy in this 
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Fig. 6. — Results of final velocity field modelling (see section 3.2.2) for the central rotating gas disk component which traces the galactic 
potential, (a - left) Best fit rotation curve (solid line) and limiting rotation curves defined by the associated parameter errors (dot-dashed 
lines) for the model with L/Lgtar = 0.25, Ropt = g pixels and inclination angle 47 ± 12°. (b - right) Enclosed dynamical mass, 

plotted as a function of radius, for the best fit model (solid lines), and within the associated parameter errors (dot-dashed lines). The 
vertical dotted lines define the region within which we can reliably trace the galactic potential, with the dashed vertical line defining the 
radius (5.25kpc) at which we determine our best-fit dynamical mass. 
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Fig. 7. — Results of our velocity field modelling (see section 3.2.2). The masked data used in the model is displayed in frame (a), the 
best fit model in frame (b), and the model subtracted residuals in frame (c). 



case.) Within a radius of 5.25 ± 1.05 kpc (i.e. 10 to 15 
pixel radius), the enclosed dynamical mass of the host 
galaxy implied by our best fit model parameters (includ- 
ing errors) is 2.05tJ;^^ x 10" Mq. 

3.3. Star formation in the host galaxy 

In addition to constraining the dynamical mass, the 
narrow Ha emission can itself be used to estimate the star 
formation rate (SFR) in the host galaxy. We measure a 
total narrow H-alpha line luminosity of 1.22 ± 0.21 x 
-[^q43 gj,gg g-i fj-onr the unmasked regions of our datacube 
(including the 12% uncertainty in the absolute flux cali- 
bration). A further maximum of 5% of the total line flux 
may have been missed in the central twenty spaxels where 
errors due to the subtraction of QSO emission dominate, 
and from the small number of other masked pixels within 
the high signal-to-noise regions of narrow line emission, 
which we add to the upper error bound. The origin of 
the Ha line emission is not necessarily solely due to the 
presence of young stars: AGN activity and shocks may 
also play an important role in the spatially extended ion- 
ization of the gas. While we cannot place our data on 
a standard BPT diagram (|Baldwin. Phillips fc Terlevich 



119811 : iVeilleux fc Osterbrock 1 119871 ) without the addi- 
tion of a second line ratio, sources with log([Nii]/Ha) 
< —0.5 are predominantly observed to be star- 
forming rather than AGN-photoioni zed sources (e.g. 
iGroves. Heckman fc Kauffmann I [20060 . Fig. |5] displays 
the unmasked regions of our datacube shaded accord- 
ing to whether the line ratio log([Nii]/Ha) is < —0.5 
or > —0.5, i.e. whether or not stellar photoionization 
is likely to be the dominant ionization mechanism, or 
whether the observed line emission is better considered 
as part of the extended narrow line region of the AGN. 
Star formation dominates in the tidal features, but is 
equally important in some of the more luminous regions 
of Ha emission towards the center of the galaxy. If the 
Ha line emission is due solely to star formation activity 
alone, the observed line flux is equivalent to a star forma- 
tion rate of 96^17 yr~^ using the scaling of lKennicuttI 
(1998). As the Ha emission can be produced by AGN 
photoionization as well as via photoionization by young 
stars, this value represents a robust upper limit on the 
overall star formation in this system. In Fig. |9K, we dis- 
play this cumulative SFR as a function of radius. In 
Fig. [9}3, we display the cumulative SFR as a function of 
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Fig. 8. — Dominant ionization mechanism pixel flagging map. 
This plot illustrates in dark grey the regions where the emission 
line ratio log([Nll]6584/Ha) < —0.5, i.e. where the emission is 
likely dominated by stellar rather than AGN photoionization. The 
light grey pixels denote the regions for which AGN photoionization 
may instead be the dominant ionization mechanism. 

radius considering only the data with line ratio values of 
log([Nii]/HQ;) < —0.5. In this case, we observe an over- 
all SFR of 50 ± lOMgyr-^ derived from Ha emission 
that is unlikely to have been photoionized by shocks or 
by the AGN, providing a useful lower limit on this quan- 
tity. The average integrated SFR per square kiloparsec 
(projected) is displayed for all data and also for the data 
with log([Nii]/Ha) < -0.5 in Fig.Et and Fig.[9li respec- 
tively. These plots clearly illustrate the importance of 
star formation throughout this galaxy, and that it is not 
restricted solely to either the inner or outer regions. 

3.4. Broad line emission and the QSO properties 

Taking the literature data for J090543. 56+043347.3, 
the mass of the central black hole has been estimated as 
9.02±1.43x IO^Mq, based on the width of its Mgii emis- 
sion l ine and the local continuum luminosity (jShen et al.l 
I2OO8I ). Additionally, th is source has a b olometric lumi- 
nosity of 3.10 X lO^^W (|Shen et al.ll2008[ ). corresponding 
to a factor of 0.25 Ledd for a black hole of this mass. 

We can also determine a black hole mass based on mea- 
surements of_the_broadHaemission line. Using the for- 
mulism of IGreene fc Hoi (f2005h . 



Mbh - (2.0ll]-j) X 10' 



L 



Ha 



10'*2ergss"\ 

2.06±0.06 



0.55±0.02 



f FWHMHa 

^ \ 103km s-1 



Mq. (2) 



We derive a spatially integrated QSO spectrum by 
combining the spectrum of the peak spaxel of our flux- 
calibrated and continuum subtracted datacube (dom- 



inated by QSO emission) with the PSF profile de- 
rived from the broad line fiux across the field of view. 
The broad Ha line has a measured FWHM of 3160 ± 
170km s~^, and an overall luminosity of 109 ± 13 x 
lO'^^ergs s~^, resulting in a black hole mass estimate of 
2.83l} i3 ^ IO^Mq. We note that the Mgii and Ha based 
black hole masses differ by a factor of ^ 3, and that the 
observed bolometric luminosity of iShen et al.l (|2008[ ) is 
equivalent to a factor of 0. 79 LEdd if we assume the Ha- 
derived black hole mass, which is likely to be subject to 
a slightly lower overall systematic uncertainty. 

4. DISCUSSION AND CONCLUSIONS 

The driving science question behind this work is to 
investigate the black hole mass vs. host galaxy mass 
relation at high redshifts, in a direct and unbiased man- 
ner as possible. Obtaining dynamical quasar host galaxy 
masses via narrow line emission kinematics is an ideal 
way forward. Although the existence of complex ve- 
locity structures can potentially present difficulties for 
such a method, our pilot object J090543.56~|-043347.3 
despite showing a non-simple velocity field can nonethe- 
less be successfully modelled. Having obtained a mea- 
sure of both black hole and host galaxy dynamical mass 
for J090543.56-h043347.3, in Fig. [11 we contrast it with 
the canonical l ow red shift Mbh vs. MsuigcDyn relation of 
iHaring fc Rixl ()2004D and also the more recent Mrh vs. 
MstRii ar and Mbh vs. MBuige,Dyn relations of Sa ni et all 
(|2011h . based on Spitzer/IRAC da ta. The dynamical 
bulge masses o f IHaring fc Rb3 (pOOl plotted on Fig. [lOli 
are either derived via first confirming that the light pro- 
files of the galaxies are bulge-dominated, and then solv- 
ing the spherical Jeans equation, or taken directly from 
the literature. This sample uses reliable pre-existing 
black hole mass estimates from the literature, derived 
from gas and stellar kinem atics and primarily sourced 
fromlTremaine et al.l (|2002[ ). For the points on Fig. [TUb . 
iSani et ain201lD carrv out 2-d image decomposition of 
the 3.6/i TO data for their sample, and then derive virial 
dynamical galaxy masses as MBuigc.Dyn = 5f7^Rc/G. 
Once again, the black hole masses for this sample are 
reliable estimates from the literature based on stellar or 
gas dynamics, or masers. Note that there is no signifi- 
cant difference between the regression lines fo r Mbh vs. 
MgjjWe.Dyn and Mbh vs. Mstdiar derived by ISani et all 
(j2011l ) from the 3.6/.t mass-to-light ratio including colour- 
correction terms. Two separate points are plotted for 
the different Mgii and Ha based black hole masses of 
J090543.56-t-043347.3. We find that, using the Mgii 
based black hole mass, the position of this object on ei- 
ther plot is a factor of ~ 2 from the mean z = relation, 
but well within the scatter. Using our Ha based black 
hole mass measurement, the position of this object lies 
directly upon the mean z = relation. This no-evolution 
result is also consistent with the theoretical modelling of 
iJahnke fc Macciol ()201lD , in which the correlation can be 
fully explained to have emerged at high redshifts from hi- 
erarchical assembly and merging statistics. Due to the 
drop both in merger rate as well as SF and BH accretion 
rate densities in the universe since z ~ 2, their explana- 
tion predicts the overall scaling relations to be largely in 
place by 2; = 1. 

It should be noted that on Fig. [10] we compare the 
dynamical mass derived for a z=1.3 galaxy with an un- 



10 



Inskip et al. 



known stellar morphology and a rotating gas disk with 
stellar masses and bulge dynamical masses in the local 
universe. Given this, it is therefore worthwhile to con- 
sider the likely subsequent evolution of this system over 
the intervening 8-9 Gyrs. The major mechanisms for 
the growth and evolution of the black hole and its host 
galaxy are BH accretion, star formation, merger activity 
and disk-to-bulge reprocessing, which we now consider in 
turn. 

The bolometric quas ar lum inosity of 
J090543.56-I-043347.3 (Shen etal] [2008h is equiva- 
lent to a black hole accretion rate of ^ 5.4Moyr~^ 
(assuming mass is converted to luminosity with an effi- 
ciency of 10%). Quasar lifetimes are typically between 
10^ a nd 10** years ( Yu & Tremainc "2002': Marconi ct alj 
[250l Hopkins & Hcrnquist 2009; KcUv ct al. 201(f7 
Taking the average of the two black hole mass mea- 
surements for J090543.56-H043347.3 (- 6 x lO^M©) 
and a timescale of 5 x 10^ years, an accretion rate of 
~ 5.4MQyr~^ would result in the central black hole 
growing by ~ 45% in total over the current active 
cycle of AGN-mode accretion. The absolute maximum 
possible growth factor for this black hole (assuming: 
[i] that it is observed at the very start of its current 
activity cycle, [ii] a maximum accretion lifetime of 10* 
years, and [iii] the lower of the two black hole mass 
estimates) is a factor of < 2, though we do stress that 
growth approaching this limiting value is unlikely. At 
the same time, it is statistically not very likely that this 
BH would undergo any further episodes of AGN activity 
between z = 1.3 and z = (gi ven the expected AGN 
fraction as a function of rcdshift: iBluck et al.ll201l[l . 

Considering star formation in the host galaxy, our ob- 
served Ha flux is consistent with a SFR of up to 100 
M0yr~^ (section 3.3). Only accounting for the line emis- 
sion which is likely due to ionization by young stars 
rather than the AGN or shocks (particularly true of the 
tidal arms and the areas of high emission line intensity 
within the host galaxy: see Fig. |3]i. Fig. ^ gives a SFR 
of 4O-6OM0yr~^. The observed star f ormation rate is sig- 
nificantly higher than average (e.g. iKarim et al.l 120111) 
for a galaxy at this redshift with a stellar mass approxi- 
mating our derived dynamical mass {2.05t.l'^4 x 10^^ Mq 
within a radius of 5.25 ± 1.05 kpc). But this may not be 
surprising given the clear presence of highly star-forming 
tidal features (Fig. |3]). 

As an aside, we note that triggering of the quasar ac- 
tivity via interactions/mergers is a possibility for this 
particular AGN, given the clearly discernable merger 
signature in the gas. The high projected star forma- 
tion rate per square kiloparsec towards the center of the 
host galaxy (within a radius of 4kpc; see Fig. [9}: and 
Fig. [9ji) might also be suggestive of links between merg- 
ers, AGN activity and circum-nuclear star formation in 
the case of this source, or of a nuclear star formation 
ring such as that observe d in the case of NGC 1097 
(Ivan de Ven fc Fathi II2O1O0 . However, the trig gering of 
AGN activity via mergers remains unproven and we are 
far from being able to e.g. determine and compare the 
age of the AG1N0 with that of stars newly formed by a 

^ The only limit we have is a mininum age of the AGN of ~ 
2 X 10^ years, from the fact that we see gas ionized by the AGN at 



merger. Overall, there is very little direct observational 
evidence that galaxy mergers trigg er AGN even in indi- 
vidua l cases at lower redshifts (e.g. iCanalizo fc StocktonI 
l2000t ). while on a population basis major merging can 
clearly be ruled out as the dominating mechanism at 
2 < 1 (|Cisternas et al.ll20lH ). 

Assuming that the local scaling relations for black hole 
mass vs. galaxy stellar mass were already approximately 
in place at z ~ 1.3, we note that mergers involving galax- 
ies with similar MeH/MTotai ratios will (statistically) 
lead to the same proportional growth for both black hole 
and host galaxy, leaving the relation unchanged. It is ex- 
pected that a stereotypical galaxy with the same mass as 
that determined for J090543. 56-1-043347. 3 would undergo 
at m ost one more major merger event (jRobaina et al.l 
120101 ). and quite likely be restricted to more minor ac- 
cretion events. Therefore, both the host galaxy and its 
central BH are likely to grow by a maximum of ap- 
proximately 50% in mass by the present day due to 
merger activity alone (s ee, e.g. iHopkins et al.l 120101 : 
Ivan Dokkum et aLll2010D . Similarly, under the assump- 
tion that the observed scaling relations for Mbh vs. Mnost 
are already in place at z — 1.3, merger activity alone 
would shift the exact placement of this object on Fig. [10] 
in a direction roughly parallel to the observed mean z = 
scaling relation. 

Even though the gas and stars trace the same gravi- 
tational potential, the gas and stellar kinematics can be 
expected to differ. While the observed gas kinematics 
for the bulk of this system are relatively settled and con- 
sistent with a rotating disk, the stellar morphology can 
not be inferred from the gas morphology, and we em- 
phasize that it remains unknown. However, given the 
observational evidence for ongoing tidal activity, even if 
this galaxy is restricted in the future solely to minor ac- 
cretion rather than major merger events, the observed 
dynamical mass at z=1.3 can be readily reprocessed into 
a stellar bulge. This object is almost certainly the pro- 
genitor of a z = bulge dominated galaxy. 

What bearing does the combination of all our obser- 
vations have on the likely subsequent evolution of this 
system and on its position relative to the local scaling 
relations? 

On the basis of the observations of 
J090543.56-I-043347.3, we can quantify the relative 
rate at which both the mass of the central black hole 
and of the host galaxy's stellar population are growing 
at the epoch at which we observe them. The annual 
increase in black hole mass of 5.4M0yr~^, as derived 
from the quasar's bolometric luminosity, is just ten times 
less than the observed (conservative) star formation 
rate of ^ 5OM0yr~^. By contrast, the scaling relations 
predict a black hole which is of the order ^ 500 times 
less massive than its host galaxy. While the central black 
hole is actively accreting, the host galaxy is forming 
new stars at a much more moderate rate relative to its 
overall stellar mass. 

Of course, we know that the duty cycle of star for- 
mation for z ~ 1.5 galaxies is far higher than the duty 
cycle of active black hole accretion. Over a single AGN 
activity cycle lasting typically 10'' to 10* years, the pro- 
portional increase in mass of the black hole is expected 

0''7 distance. 
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Fig. 9. — (a - top left) Cumulative SFR (solid line) plus error bounds (dashed lines) as a function of radius, calculated using the scaling 
of IKennicutH l|1998l ). (b - top right) As frame (a), but only including data from spaxels with an emission line ratio log([Nll]6584/Ha) 
< —0.5, i.e. where the emission is likely dominated by stellar rather than AGN photoionization (see fig 8.). (c — lower left) Average SFR 
per square kiloparsec (projected) as a function of radius (solid lines), plus error bounds (dashed lines), (d - lower right) As frame (c), but 
only including data with log([Nll]6584/Ha) < -0.5. 



to be ^ 45%, as noted above. A similar proportional 
increase in the stellar mass of the host galaxy can be 
achieved if the observed level of star formation persists 
over a substantially longer timescale than the AGN activ- 
ity: a galaxy of 2 x IO^^Mq with a star formation rate of 
~ 5OM0yr~^ would increase in stellar mass by 45% over 
the course of ~ 2Gyr. In such circumstances, the black 
hole mass vs. host galaxy mass ratio of this object would 
remain consistent with the local relation. However, we 
note that the required star formation timescale of several 
Gyr is well beyond the typical merger timescales. 

An alternate approach to constraining the expected 
future increase in stellar mass of J090543. 56+043347. 3 
with time until z = is to consider star formation rates 
for the galaxy population as a w hole. More speci fically, 
following the cosmic average (e.g. iKarim et al ."2011) the 
expected levels of star formation activity would result in 
a growth of typically a factor of about 1.5 to 3 in stellar 
mass between z = 1.3 and the present day, depending on 
the exact details of the subsequent star formation activ- 



ity within the galaxy. 

We can combine this empirical estimate of the stellar 
mass growth through star formation activity with the 
increase in mas s due purely to the e xpected merger ac- 
tivity (~ 50%; iHopkins et all 120101 ). OveraU, the host 
galaxy is very likely to more than double in mass (in- 
creasing by > 0.7 dex in log(Mtotai)) between z — 1.3 
and z = 0. As the majority of any disk stellar mass 
will have been re-processed into the galactic bulge via 
galaxy mergers or tidal interactions, we would expect 

log(Mtotal,z=o) « log(MBulgG,z=o)- 

Similarly, we can combine the expected growth in BH 
mass due to AGN accretion and merger activity to pre- 
dict the overall increase in Mbh by z = 0. These 
processes are likely to result in an overall increase of 
< 0.25 dex in log(MBH) for a system such as this; taking 
our less-likely absolute upper limit for BH accretion, the 
maximum increase in Mbh is ^ 0.7 dex. 

Any change in position of this object by z = rel- 
ative to the local relation will therefore most likely be 
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Fig. 10.— J090543. 56+043347.3 and the black hole mass vs. host 
galaxy mass relation (a — top) Position of J090543. 56+043347. 3 
(stars) relative to the z = Mbh vs. MBuigo,Dyn correlation, with 
datapoints from lHarine fc Hixl 120041) . The upper (blue) star uses 
the black hole mass derived from the Mgll emisison line, while the 
lower (red) star uses our Ha based black hole mass. Note that the 
total dynamical mass is plotted for J090543. 56+043347. 3, and that 
this represents an upper limit to the bulge mass at z = 1.3. The 
solid line displays the bisector linear regression fit of Haring & Ri:)3 
((2004), while the dotted line gives the linear regression fit for the 
Mbh vs. Mstcllar relation derived bv Sani et al t201l] ). (b - bot- 
tom) Position of J090543. 56+043347. 3 (stars, color-coded as in 



frame (a)) relative to the 2 = Mbh vs. M 



Bulge, Dyn 



correla- 



tion, with datapoints and linear regression fit (dashed line - note 
that this is indistinguishable from the Mbh vs. MsteHar regressio n 
fit denoted by the dotted line on frame (a)) from lSani et afll2011l ). 



dominated by changes in the host galaxy mass. Overall, 
regardless of the exact details of the further evolution 
of this object, it is most likely that the average of the 
two points plotted for this object on Fig. [10] will move 
to a position which is even closer to the best fit z=0 re- 
lations displayed in Fig. IIOI and in any case will remain 
consistent with the observed scatter. 

In conclusion, we return to the observational aspects 
of the work presented here. We believe that we have 
presented a study that is ground-breaking in the respect 
that it demonstrated how AO-assisted near-IR integral 
field spectroscopy can be successfully applied, and used 
to determine the dynamical mass profile for the host 
galaxy of a typical quasar. The observed BH masses and 
overall dynamical mass of the host galaxy of this object 
at z = 1.3 are consistent with the observed present-day 
relations for Mbh vs. MBuigo,Dyn and Mbh vs. Mstdiar, 
and aside from a potential conversion of disk- to bulge- 
mass through minor mergers or the likely ongoing ma- 
jor merger, no significant evolution in the black hole vs. 
stellar mass relationships between z = 1.3 and z = 
is required or expected. Provided that one takes into 
account the minimum extended emission line fiux re- 
quirements and the likelihood of favourable gas kinemat- 
ics and orientation in the observed velocity field, this 
particular diagnostic route is widely applicable to other 
quasars in the redshift range 1 < z < 3, at the peak 
epoch of quasar activity. Analysis of a further seven 
objects in our pilot sample is currently underway. Al- 
though J090543. 56+043347.3 itself displays a relatively 
impressive star formation rate, this diagnostic method is 
clearly feasible at much lower narrow Ha fiuxes, while al- 
ternatives such as CO-based dynamical masses are most 
successful for extreme, atypical objects. 
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